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It is widely believed that the critical behaviour of the two-dimensional Ising model 
is only slightly modified by the introduction of non-frustrated disorder [|l], |2j. Such 
changes are given by logarithmic corrections to pure-system power-law singularities. In 
terms of a field-theoretical (or renormalisation-group) description, disorder is said to be a 
marginally irrelevant operator M. Specific forms have been proposed, and numerically 
tested, for the corrections to bulk quantities as specific heat, magnetisation and initial 
susceptibility [[I]]. The overall picture emerging from such analyses tends to confirm 
predictions of log-corrected pure- Ising behaviour. Early proposals of drastic alterations 
in the values of critical indices || have thus been essentially discarded. However, recent 
results [§, have appeared, according to which critical indices would vary with disorder, 
but so as to keep the ratio 7/1/ constant at the pure system's value (the so-called weak 
universality scenario ||). In order to solve this controversy, it is important to undertake 
independent tests of several aspects of the problem. In the present work we examine the 
finite-width corrections to the first gap of the Lyapunov spectrum of the disordered Ising 
model on strips. This gap is related to the typical (as opposed to averaged) behaviour 
of spin-spin correlation functions, as explained below. 

While the average of a random quantity Q is simply its arithmetic average over 
independent realisations Qi, Q = (1/N) Y^iLi Qi, its typical (in the sense of most 
probable) value is expected to agree with the geometrical average: Q typ = Q m . P . = 
exp (in Q] || [10], [TTJ |I2|| . Depending on the underlying probability distribution, 
Q and Qtyp may differ appreciably, as is the case when one considers correlation 
functions [g [TJ. 

It has been predicted, on the basis of field-theoretical arguments ||, that in a bulk 
system (as opposed to the strip geometry used here) the typical decay of spin-spin 
correlation functions at criticality on a fixed sample is given by 

(Widtyp oc iT 1/4 (A In i?)- 1/8 (fixed sample) (1) 

for ln(Alni?) large, where A is proportional to the intensity of disorder. This way, 
pure-system behaviour (power- law decay against distance, (cxocxr) ~ i? _r? with 77 = 1/4) 
acquires a logarithmic correction. On the other hand, if one considers the average, over 
many samples, of the correlation function, logarithmic corrections are expected to be 
washed away || resulting in a simple power-law dependence: 

Wo^r) ^ -R -1 ^ 4 (average over samples). (2) 

The distinction between typical and average correlation decay was not explicitly 
discussed in the early field-theoretical treatment || which predicted (cr cr^) oc e - A ( l]0 -R) 2 ^ 
a sort of behaviour for which no evidence has been found in subsequent investigations 
Recent numerical work claiming weak universality to hold || M does not address the 
issue either, though it is easy to see that the procedures used in those calculations 
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pick out averaged correlations, as they rely respectively on variants of the fluctuation- 
dissipation theorem H] or on explicit averaging 0. 

In Ref. J7|, an analysis of phenomenological renormalisation estimates of the 
correlation-length exponent v seems to point against logarithmic corrections to its pure- 
system value, and in favour of a disorder-dependent exponent. However, that is based 
on data for very narrow strips [L < 8) and mostly relies on trends apparently followed 
at weak disorder. Though one series of moderately-strong disorder results is exhibited as 
well, no attempt is made to fit either set to the form predicted by theory, which crucially 
includes a disorder-dependent crossover length Q. Indeed, a systematic treatment of 
the averaged correlation lengths at the exact critical point, spanning a wide range of 
disorder and watching for the disorder-dependent crossover mentioned above, eventually 
uncovers the expected logarithmic terms ||. Turning back to the exponent rj, recall that, 
as far as dominant behaviour is concerned, both weak- and strong-universality concur 
in predicting rj = 1/4. Thus, our strategy here will be to analyse the subdominant 
terms. Our goal is to show that the finite-size corrections to the typical ratio of decay 
of correlations behave consistently with what is expected when a marginal operator is 
present (see below). 

In numerical simulations, one considers finite lattices (L x L) or finite-width strips 
(LxN, N —>■ oo) and sets the temperature at the critical point of the corresponding two- 
dimensional system. For the nearest-neighbour random-bond Ising model on a square 
lattice, with a binary distribution of ferromagnetic interactions 

P(^) = \(S(J lJ -J ) + S(J lJ -rJ )) , 0<r<l , (3) 

the critical temperature (3 C = l/ksT c is exactly known [13|, [14| from self-duality as a 
function of r through: 

sinh(2& J ) sinh(2/5 c r J ) = 1 . (4) 



From Monte Carlo work on L x L random-bond systems ||15| , it has been found that 
the average correlation function at criticality is numerically very close to the exactly- 
known |16| value for a pure system of the same size at its own critical point, thus 



providing evidence in favour of Equation (Q). Similar conclusions have been drawn 



for the corresponding quantities evaluated on strips [17], |18| , where the exact critical 
correlation functions for the pure Ising model are known from conformal invariance [19[ . 

Here we provide a test of the consequences of Equation (Jl]), when correlations are 
considered on a strip geometry. In this case, contrary to that of Equation (|2]), no exact 
finite-L results are available for comparison; thus one must resort to finite-size scaling 
concepts |]2Tfl , in order to unravel signs of the expected bulk behaviour from trends 
followed by finite-system data as L — > oo. 

The procedure used is as follows. It is known that the typical, or most probable, 
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spin-spin correlation function on a strip decays as 

(cT CrR)typ OC exp (-R/Ztyp) , £typ = A L ~ A L , (5) 

where A° L and are the two largest Lyapunov exponents of the (random) transfer 
matrix for a strip of width L |TU| , |Tl"| , |T2[. On the other hand, conformal invariance 
predicts [|21| that, when the Hamiltonian of a homogeneous two-dimensional system 
contains a marginal operator, the spectrum of eigenvalues E n of the transfer matrix on 
a strip is such that 

E n -E = (2tt/L) (x n + dj In L) + . . . , (6) 

where x n is the corresponding scaling dimension, d n is an n-dependent constant and 
periodic boundary conditions are used across the strip. Since (i) disorder in the 
two-dimensional Ising model is believed to be a marginal operator, (ii) Lyapunov 
exponents of transfer matrices in random systems are the counterparts of eigenvalues in 
homogeneous ones, and (Hi) numerical evidence shows that conformal-invariance results 
derived for uniform systems may be extended to random cases, provided that suitable 
averages are taken 0, |17] , we shall examine sequences of estimates A° L — for varying 
L and try to fit them to Equation ((^) with n = 1 and 2x\ = rj = 1/4 fl9|| . 

We have used strips of width L < 13 sites, and length N = 10 5 for L = 2 — 11 
and 5 x 10 4 for L = 12 and 13. Two values were taken for the disorder parameter r 
of Equation ([3]): r = 0.5 and 0.01. Details of the calculation are given in Ref. fT7 



where the data used here are displayed as well, in plots of (L/n)(A° L — A^) against 
1/L 2 . The latter variable was used because it is exactly known [E2J that, for strips of 
pure Ising spins, the leading corrections to the L _1 dependence of £ -1 are proportional 
to L~ 2 . Our purpose there was to show that the correlation length £ ave coming from 
direct average of correlation functions, (o"o<7r) ~ exp(— R/^ ave ), indeed scales as its 
pure-system counterpart, while £ typ does not. At the time we did not investigate the 
behaviour of £ typ in detail, though it was noticed that for strong disorder the curvature 
of plots of L/n^typ against only became smaller for close to zero. 

In Figures 1(a) and 1(b) we show L/ir£ typ = (L/7r)(A° — A^) against 1/lnL, 
respectively for r = 0.5 and 0.01. For weak disorder r = 0.5 the fit to Equation (|6|) 
is very good, with 2x i = rj = 1/4. While for r = 0.01 strong fluctuations increase 
the amount of scatter, the overall picture still is consistent with Equation (|6|). Table 
1 shows results from least-squares fits of data in the range L = 6 — 13. This interval 
was chosen in order to minimise the accumulated standard deviation \ 2 P er degree of 
freedom [|17]]. The extrapolated r] is expected to be universal as long as r ^ 0; our error 
bar for r = 0.01 indeed includes 7] = 1/4, though admittedly at the edge. On the other 
hand, d\ clearly changes, increasing with disorder. The situation differs from that of 
the g-state Potts model with q = 4 [ZT[ ^3], |24| where there is no continuously tunable 
parameter, as marginality depends on q — 4 being strictly zero. 
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Table 1. Results from least-squares fits to Equation (jfj), for L = 6 — 13. 





r 


V 




x 2 


0.5 
0.01 

Expected 


0.250 ± 0.006 
0.29 ± 0.04 

1/4 


0.018 ± 0.013 
0.38 ± 0.09 


.021 
.17 



As an additional check on the ideas developed above, we examined data for the 
probability distribution of critical spin-spin correlation functions [pl| for r = 0.25. 



For spin-spin distances R = 5 and 20, and strip widths L = 3, ... 13 we picked the 



averages \nG(R) = In ((coOr)). The quantity explnG(i?) is expected [JTTJ, [12] to 
scale as (cr aR) t yp- For each L the slope of a two-point semi-log plot of explnG(i?) 
gave an approximate value for From a plot of L/7i^ typ against 1/lnL one gets 

r] ~ 0.26, di ~ 0.03, consistent with the above values derived directly from Lyapunov 
exponents, and with the assumption that d\ varies continuously (and monotonically) 
against disorder. 

We have analysed numerical estimates for the first gap of the Lyapunov spectrum 
of the self-dual random-bond Ising model on strips. We have shown that finite-width 
corrections can be fitted very well by an inverse logarithmic form, predicted to hold 
when the Hamiltonian contains a marginal operator. The present results contribute, 
albeit indirectly, to the growing body of evidence favouring strong universality (that 
is, pure-system behaviour with logarithmic corrections) in the two-dimensional random 
Ising model fl], §, p|, £§ . This is to be contrasted to recent work || 0] , according to which 
critical indices would vary with disorder, but so as to keep the ratio 7/V constant at the 
pure system's value (the so-called weak universality scenario |§). 
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Figure captions 



Figure 1. 77 = L/n^typ, with £ t ~p = — of Equation (g), against 1/lnL. The 
square on the vertical axis is at the pure system value rj = 1/4 . Straight lines are 
least-square fits for L = 6 — 13. (a) : 7- = 0.5; (b) : r = 0.01 . 
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